PERIODIC BILLIARD ORBITS ON n-DIMENSIONAL 
ELLIPSOIDS WITH IMPACTS ON CONFOCAL 
QUADRICS AND ISOPERIODIC DEFORMATIONS. 
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Abstract. In our paper we study periodic geodesic motion on 
multidimensional ellipsoids with elastic impacts along confocal quadrics. 
We show that the method of isoperiodic deformation is applicable. 



1. Introduction 

One of the classical integrable problems is the geodesic motion on 
the ellipsoid. Complete integrability of the 3-dimensional problem was 
proven by Jacobi [25] • Explicit formulas in terms of genus 2 6'-functions 
were found by Weierstrass [H]. Integrability of the multidimensional 
problem was established by Moser [32]. Connections with the Neumann 
problem were exploited by Knorrer [27]. Integrability of the Neumann 
problem has a very transparent explanation, based on the following 
observation by Trubowitz and Moser (see [32]): taking properly nor- 
malized eigenf unctions at the right ends of gaps for the (A^ — l)-gap 
1-dimensional stationary Schrodinger operator we obtain the coordi- 
nates of the particles. 

The isomorphism between the geodesies on the ellipsoid and the Neu- 
mann problem is nontrivial, because it includes a nonlinear change of 
independent variable. The motion on the Jacobian for the Neumann 
problem is uniform, and the motion for the geodesies on the ellipsoid 
is along straight lines, but the speed in the natural parameter is not 
constant. It means, that the analog of the Trubowitz-Moser isomor- 
phism maps this dynamics into the Harry-Dym equation instead of 
the Korteweg-de Vries (KdV) equation. On the level of the associated 
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spectral problems it means, that one has to consider the 1-dimensional 
acoustic problem 

(1) -d'f^{E,t) = Eu{t)^{E,t), E = l/X 
instead of the 1-dimensional Schrodinger problem 

(2) - dl^{E, a) + u{a)^{E, a) = E^{E, a). 

The connection between the the geodesies on the ellipsoid and the 
Harry-Dym equation was first observed by Cao |8], and independently 
by Veselov [H], [H. 

From the geodesic motion on the ellipsoid one can construct an inte- 
grable billiard motion taking the geodesic dynamics on the n-dimensional 
ellipsoid 

(3) Q=\^ + ... + ^ = 

I Oi CLn+l 

combined with elastic reflections along one or more confocal quadrics 

Q.. = (^V + --- + ^^ = 4' ^ = i'---'-- 

dg being arbitrary fixed parameters. An amazing property of this mo- 
tion is that all conservation laws remain invariant after reflections, 
therefore this dynamical system is also completely integrable. Such 
property was first observed by Chang and Shi [9] for the triaxial ellip- 
soids, while its multidimensional generalization was obtained in |lUj. 
As a consequence the trajectories between impacts (arches of geodesic) 
or their continuations are tangent to the same set of — 1 quadrics 
(caustics) in R""*"^ that are confocal to Q and Qds- the same papers 
it was proven the Poncelet property: if one trajectory is closed after p 
bounces, all the trajectories sharing the same constants of motion are 
also closed after p bounces and have the same length. 

The billiard on the elhpsoid is described by the maps Eg : (x*^^^ , v^^^}) y 
(x(^+^),wi^+^^), where x^^\v^Z} E R"+\ N e N, are respectively 
the coordinates of the A^-th impact point on Q fl Qd^ and the out- 
going velocity at this point, whereas x''^~^'^\vl^^^^^ denote the same 
objects at the next impact point. This billiard map is transcendental 
in the following sense: the shift takes place on a non-Abelian variety 
and, as a result, it is not fixed. This property was first observed by 
Veselov [IQ] for the billiard on triaxial ellipsoids and generalized for 
higher dimensions by Abenda and Fedorov in [3J. More precisely, the 
generic complex invariant manifolds are not Abelian varieties, but open 
subsets of theta-divisors of n-dimensional hyperelliptic Jacobians of a 
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given hyperelliptic curve T. Thus, in the real domain, the new values 
(^^(Af+i)^ g^j,g determined by solving a system of transcendental 

equations which involve the inversion of hyperelliptic integrals, whereas 
in the complex coordinates the billiard map is infinitely many valued. 

We stress that all of these properties stand in sharp contrast to 
both the Birkhoff (p = 0) [3 [Ml EH El] and the Hooke (p = 2) [H] 
billiard systems, where the billiard map is given by a shift by a constant 
vector on a complex invariant manifold which is an open subset of 
an Abelian variety. Due to the above properties, the billiard maps 
Bs can be regarded as a discrete analog of the algebraic integrable 
systems with deficiency studied, in particular, in [l3l [21 [H [H]). The 
latter systems, although being Liouville integrable, are not algebraic 
completely integrable: their generic complex solutions have movable 
algebraic branch points and are single-valued on an infinitely sheeted 
ramified covering of the complex plane t, so they possess the so called 
weak Painleve property. The branching of the billiard maps Bg can be 
viewed as a discrete version of such a property (see also [2T]). 

The development of the finite-gap approach to the 1-d Schrodinger 
operator and the KdV equation was started by Novikov ^34j- Generic 
finite-gap potentials are quasiperiodic. Selection of purely periodic po- 
tentials is rather nontrivial. 

The principal periodicity constraints discussed in the literature are 
the following: 

(1) Periodicity in one real variable. 

(2) Periodicity in two real variables. 

(3) Double periodicity in one complex variable. Such solutions are 
called elliptic solitons. 

The problem of selecting space-periodic potentials for the 1-D Schrodinger 
operator was solved by Marchenko and Ostrovski |3T] in terms of con- 
formal maps. It is also possible to study the variety of all spec- 
tral curves generating purely periodic potentials, using the so-called 
isoperiodic deformations (they can be interpreted as Loewner equa- 
tions for these conformal maps). Period preserving deformations first 
arose in the papers by Ercolani, Forest, McLaughlin and Sinha [T7], 
Krichever [3Dj, Schmidt [35j. Grinevich and Schmidt ^22j suggested 
to use isoperiodic deformations for selecting pure periodic solutions, 
they also showed that these flows can be transformed to a very simple 
algebraic form by extending the phase space. 

The finite-gap theory for the acoustic problem ([1]) was developed by 
Dmitrieva and collaborators [11], [12], [T3]. This problem can be trans- 
formed to the Schrodinger problem ([2]) by a reciprocal transformation 
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over the independent variable. Equations for the isoperiodic deforma- 
tions of this problem were written by D.Zakharov [IS]. These results 
can be automatically applied to the periodic geodesies. 

The first non-trivial examples of elliptic solitons for the KdV equa- 
tion (and 1-D Schrodinger operators potentials respectively) were found 
by Dubrovin and Novikov [15]. These solutions correspond to genus 
g = 2 curves. Airault, McKean and Moser [6] established the con- 
nection of KdV elliptic solitons with a special reduction of the elliptic 
Calogero-Moser system valid for arbitrary genera. An analogous con- 
nection between the Kadomtsev-Petviashvili (KP) equation and the 
full elliptic Calogero-Moser system was discovered by Krichever [29]. 
A complete characterization of the double-periodic in the complex x- 
variable KdV and KP solutions in terms of so-called tangential cov- 
erings was obtained by Treibich and Verdier [36j, t37j. Due to the 
Trubowitz- Moser isomorphism these results can be also applied to fully 
characterize and construct real geodesies and real billiard trajectories, 
which are double-periodic in the complex length [TH], [3], [I]- 

It is natural to pose the problem of constructing real billiard tra- 
jectories on ra-dimensional ellipsoids which are simply-periodic. In our 
paper we show, that the isoperiodic technique can be applied to this 
problem. 



In this section, to settle notations, we briefly outline the description 
of the geodesies and billiard motions on a quadric. The geodesic motion 
on an n-dimensional ellipsoid Q is well known to be integrable and to 
be linearized on a covering of the Jacobian of a genus n hyperelliptic 
curve [32]. Namely, let t be the natural parameter of the geodesic and 
Ai, . . . , A„ be the ellipsoidal coordinates on Q defined by the formulas 



2. Geodesics on Ellipsoids 



(4) 




I = 



n + 1. 



Then, denoting A^ = d\k/dt the corresponding velocities, the total 
energy - (x, x) takes the Stackel form 



n 
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According to the Stackel theorem, the system is Liouville integrable. 
After fixing the constant of motion c^, = 1, . . . , n — 1 and assuming 
if = 1/2 we obtain Dubrovin equations: 

(5) A. = -2^^ . 

Afc n (Afc — Aj) 

Here 
(6) 

wl = -Afc(Afc-ai) ■ ■ ■ (Afc-a„+i)(Afc-ci) ■ ■ ■ (Afc-c„_i), = 1, . . . ,n. 

The evolution of the divisor A^ is defined on the genus n hypereUiptic 
curve 

71+1 71—1 2?! 

(7) r : {w' = -Xl[{X-a,)l[{X-c,)} = {w' = -l[{X-b,)}, 

i=l k=l i=0 

where we set the following notation throughout the paper 
(8) 

{ 0, ai < ■ ■ • < a„+i, ci < ■ ■ • < Cn-i } = {bo = <bi < ■■■ < fean}- 

The constants of motion have the following geometrical mean- 
ing (see [32]): the corresponding geodesies are tangent to the quadrics 
Qci, • • • , Qc„-i of the confocal family Qc = {xl/{ai - c) H h x^_^i/(a„+i - c) = l}. 

In particular, the reality condition for geodesies on ellipsoids is equiv- 
alent to either q = b2i or q = 621+1, i = I, . . . ,n — 1 (see [261 [5]). Since 
we are interested in the reality problem, we take,without loss of gener- 
ality, b2i-i < Xi < b2i, i = 1,. . . ,n. 

Let us recall (see [8], [H]) how to reconstruct the geodesic motion 
from the solutions of the acoustic problem ([1]). The n-gap eigenfunc- 
tions of ([1]) are defined by the following properties: 

(1) For any t, '${'y,t) is meromorphic in 7 on the spectral curve 

r\{A = o}. 

(2) For any t, \1'(7, t) has n simple poles at the points 71(0),. . . ,7ri,(0). 
Here Xj{t) is the projection of 7j(t) to the A-plane. 

(3) The there exists a function cr(t) such, that 

is regular near the point A = 0. 

(4) 

(9) ^(A,t) = 1 + ^ + 0^ as A ^00. 
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Let US define 



(10) 



Xj{t) = ^(aj,t) 



% n - At(o)) 

fc=i 

n - (^k) 



Consider the following pair of meromorpliic forms in the A-plane: 



n(A-A,(0)) 



:ii) 



fii = ^(7,t)^'(r7,t) 



fc=i 



n+l 

n (A - a,) 

k=l 



-d\. 



'12) 



^2 = ^(7,t)^(r7,t) 



An(A-Afc(0)) 

k=l 

n+l 

n (A - a,) 

k=l 



dX. 



Here r denotes the hyperelliptic involution. From it follows, that 
the integrals of Qi, Q2 over sufficiently large contour are equal to 27ri. 

The residues of fii and ^2 at the point ^— and ij(t) respectively. 

Therefore we obtain ([3]) together with the normalization condition: 

(13) E*?W = 1- 

The motion is governed by the Lagrange multipliers principle: 



Xj{t) = —u{t)- 



After the re-parameterization 



Ai ■ ■ • A„ 



: da. 



(14) dt ^ 

Vbib2 ■■■hn 

(this re-parameterization coincides with the change of independent 
variables, connecting the acoustic operator ([1]) with Schrodinger op- 
erator ([2])) the evolution of Aj is described by quadratures 

(15) 

r\fc-i 



29 

The quadratures involve n independent holomorphic differentials on 



-^dXi + ■■■ + -^dXn 

2Wi 2Wn 



da for k = 1 , 
for A; = 2 n, 



r. 



PERIODIC BILLIARD ORBITS ON n-DIMENSIONAL ELLIPSOIDS 



7 



Let ai, i = 1, . . . ,n he the conventional homological basis depicted 
in Figure 1 and 



(16) u 



\ TT^i ' ^ 



n, 



be the usual basis of holomorphic differentials. ( ITSll give rise to the 
Abel-Jacobi map of the n-th symmetric product r^") to the Jacobian 
variety of F (see for example |23j): 

Pi Pn ( 

r r \ 2a + const., for j = 1 , 

^^^^ J '^^ ^ ^ J I const., for j = 2, . . . , n, 

Po Po ^ 

where and Pq is a fixed basepoint and Pk = {Xk, Wk) G F, = 1, . . . , n. 
Then, the geodesic motion in the new parameterization is linearized on 
the Jacobian variety of F. Its complete theta-functional solution was 
presented in [l^ (n = 2), and in [26] (n > 2), whereas a topological 
classification of real geodesies on quadrics was made in [5] and an alge- 
braic geometric characterization of a dense set of real closed geodesies 
on ellipsoids has been given in pQ. 




The billiard map is associated to the geodesic motion of a point on Q 
with elastic bounces along a set of confocal quadrics Qd^, s = 1, . . . ,r, 
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di < d2 < . . . < dr. We are interested in real motions only, and we 
assume, that all reflectors Q fl Qd, are met at least once. Therefore we 
have the following constraints to the positions of the reflectors: 

n 

(1) The points ds are located inside gaps: dg E [J ]&2j-i, &2j[, s = 

i=i 

l,...,r 

(2) Each gap [&2j-i,&2j] contains at most 2 points ds- 

(3) If the gap [b2j-i,b2j] contains exactly 2 points dg, dg^i, then 
the projection of the divisor point Xj{t) is located inside the 
interval [ds, ds+i]- 

(4) If the gap [62j-i,&2j] contains exactly 1 point dg, then the pro- 
jection of the divisor point Aj(t) remains forever inside one of 
the following intervals: [62^-1, c^s] or [ds,b2j]- 

The divisor dynamics can be described in the following way. It is 
guided by the Dubrovin equations (jS]) between the subsequent collisions 
with the reflectors. If a divisor point Xj collides with reflector dg, it 
"jumps" to the opposite sheet of T with respect to the hyperelliptic 
involution: 

(18) T : {\j,Wj) ^ {\j,-Wj). 

Each point dg has two preimages in the spectral curve F: = 
{dg, +w{dg)), D~ = {dg, —w{ds)), w{dg) > 0. Let us introduce the 
following notation: 



(19) D': 



(20) 

D[ 

Prom the Dubrovin equations (Q it follows, that the divisor points 
always jump from D'^ to 

Under the Abel-Jacobi map ( fTTl) . the condition Xj = dg defines two 
codimension one subvarieties and G°"* in Jac(r), and each reflection 
on the Q fl Qd^ corresponds to the jump from one of these subvarieties 
to the other one. 

Hence, for fixed constants of motion, the coordinates x of impact 
points on Q n Qd^ with velocities v are described by degree n — 1 
divisors {Pi, . . . , Pn-i}, that is, by a point ip on 9™ or 6°*^*. 



Dt if 



(Aj - dg) 
n (A, - Xk) 



> 



D- if J^l^^<o 

n (A, - A.) 
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Let US denote 



(21) 



n\T 



Din 



, UJ. 



Let us recall the situation in case of one reflector di = d, = 



0? 



out 



0c 



D\^^^ q = qi. The bilhard map is 



transcendental, as it was flrst observed by Veselov for n = 2 [40j. An 
algebraic geometric description of the motion between impacts and at 
the elastic bounce along Q H is given by the following proposition 

Theorem 3.1. The geodesic motion on Q between subsequent impact 
points with coordinate x^^^ and x^^^^-* on Q D Qd corresponds to the 
straight line uniform motion on Jac(r) between 0°"* and along the 
ui-direction. The point of intersection with 0~ gives the next coordi- 
nate x^^^^^ and the ingoing velocity v[^^^\ At the point x'^^^^^ the 
reflection vl^^^"^ ^ ^out'*'^^ (from the ingoing to the outgoing velocity) 
results in jumping back from to 0°"* by the shift vector q, which 
does not change x^^^^-* . Then the procedure iterates. 

The process is sketched in Figure 2. Here A^^ denotes the image 



through the Abel map of the divisor, associated to the point (x 



{N) AN) 



Aq^^. denotes the image of (x 



^out 



respectively: 



.1 




out 



out 



Ur 



Figure 2: The billiard dynamics in the Jacobian. 
To study the billiard dynamics it is important to choose an appropri- 
ate canonical basis of cycles. We assume that F has the following sys- 
tem of cuts: an inflnite cut ] — oo, 0] and flnite cuts over gaps [&2j-i, b2j] 
(all cuts are real). We choose the a-cycles in the standard way (see 
Figure 3). It is important that they do not intersect the cuts, and 
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consequently, the trajectories of the divisor. If the gap [62^-1, &2j] does 
not contain reflectors, the corresponding /3j-cycle is also chosen in the 
standard way, and it intersects exactly once the cut [&2j-i5&2j]- If the 
gap [b2j-i,b2j] contains at least one reflector, we choose the cycle Pj 
so, that it does not intersect the divisor trajectory. In the case of two 
reflectors in the gap [627-1, &2j] there are two possible choices of (see 
cycle /?3 in Figure 3). 




Figure 3: The canonical basis of cycles in the presence of reflectors. 

Let us us also introduce an auxiliary set of curves ai, . . . , a„. If 
the gap [b2j-i,b2j] does not contain reflectors, aj is the preimage of 
[b2j-i,b2j] with the proper orientation and homologically aj ~ aj. If 
the gap [&2j-i, b2j] contains reflectors, aj is the part of the preimage of 
[b2j-i,b2j] where the divisor moves. The curves cxj inherit the orienta- 
tions from the basic cycles aj. 

Let us define the auxiliary basis ui,. . . , a}„ of holomorphic differen- 
tials with the following non-standard normalization: 

(22) j '^''^ 

Proposition 3.2. A basis Ui,. . . , cj„ with normalization [2^) exists 
and is unique. 

The existence and the uniqueness of such basis of holomorphic dif- 
ferentials, follows immediately from the following lemma. 

Lemma 3.3. Let Q^i^, . . . , Q(n), G = be a set of paths in the 

cycles ak such, that 
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(1) The velocity vector for all paths ^(j) is everywhere non-zero. 

(2) All paths Q(j) belonging to the same have the same orien- 
tation with respect to the standard orientation on a^. Let us 
denote: Gk = G ^ oik, i-^- Qk is the union of all trajectories, 
lying in the cycle a^. 

Then for any k such, that Qk ^ there exists a holomorphic differential 
uj such, that 



(23) iu = 5. 



>kl- 



Qi 



Proof of Lemma I3.3t Let us consider the following holomorphic dif- 
ferential 

n(A-r/,)rfA 

(24) u = C'-^ , 

w 

where fjj G aj. Differential u has a constant sign on ak, therefore 

(25) juj^O. 

If g. = then 



(26) Ij = J ^ = ^- 

Gj 

If Qj 7^ and fjj lies in the end of the interval [&2j-i, &2j] then integral 
Ij is nonzero. The signs of Ij for fij = &2i-i and fjj = b2j, are opposite, 
and they do not depend on the position of other zeroes fji. Therefore 
we can tune the set fjj so, that all Ij = for all j ^ k. To complete 
the proof is it sufficient to choose a proper normalization constant C. 

Remark 3.4. This proof uses the following topological Lemma: Let f = 
(/^(xi, . . . , Xn), . . . , /"(xi, . . . , Xn)) be a map from the n-dimensional 
cube J"; —1 < Xk < 1, k = 1, . . . ,n to the n-dimensional space M" 
with the following property: for all k f^{x) < at the face Xk = —1 
and f^{x) > at the face Xk = 1- Then the preimage of the origin 
X = (0, . . . , 0) is non-empty. 

The proof of the Lemma is based on standard topological arguments 
and we do not present it here. 
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Corollary 3.5. Let Coi,. . . , Un be the standard basis of holomorphic 
differentials 



(27) / uok 



-'jk, 



J be a subset of the set {1,2, ... , n} such, that for any k E J Qk ^ ■ 
Then the matrix 

(28) Aki = [ooi, k,lej 



is non-degenerate. 

To prove this Corollary it is sufficient to assume, that for all k such, 
that I < k < n, k ^ J' Qk = ctk and to apply Lemma 13.31 The 
differential obtained has zero integral over all contours ak, k ^ J , 
therefore it is a linear combination oi d^i, I & J'. 

Periodic billiard trajectories. First of all, let us recall, that there are 
two different descriptions of periodicity for the billiard motion. If the 
billiard trajectory is closed, it implies, that the divisor motion is peri- 
odic with the same period. From the periodicity of the divisor motion 
it follows, that the motion on the ellipsoid is also periodic, but the 
period may be greater, since the same divisor corresponds to 2"+^ dif- 
ferent points in the phase space assuming all constants of motion to be 
fixed. These points can be obtained one from another by the reflections 
Xj —Xj. We see, that each periodic divisor motion corresponds to 
2n+i cliffg];gnt periodic trajectories. 

In the following we study necessary and sufficient conditions for the 
periodicity of the divisor motion. 

Let us recall the necessary and sufficient condition on the spectral 
data to generate periodic billiard trajectories. In our setting for arbi- 
trary number of reflectors we have: 

Lemma 3.6. The billiard trajectory is periodic if and only if there exist 
a collection of integers mi,. . . ,mn such that 

n „ 

(29) ^rrij ujk = 0, k = 2,...,n. 

In such a case the length of the billiard trajectory in the rescaled variable 
a is given by 

n „ 

(30) ^ nij uji = 2L. 

i=i - 
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Here Uk is the basis of holomorphic differentials ( |T6i) normalized by 
the asymptotics near the Weierstrass point A = 0. 

The above conditions were also given in [H] in an equivalent for- 
mulation and are in general transcendental in the parameters of the 
problem (the square semiaxes, the constants of the motion and the po- 
sition of reflectors). In the case of one reflector, in [3] the authors gave 
a sufficient condition so that the periodicity condition for the billiard 
trajectory is algebraic in the parameters. 

To apply the method of isoperiodic deformations we need an alter- 
native characterization of the periodic spectral data. The remaining 
part of this Section is dedicated to this alternative characterization in 
terms of the so-called generalized quasimomentum dp. 

First of all we assume that our billiard motion is periodic, and we 
construct the generalized quasimomentum, satisfying the properties 
1-9 in Lemma I3.10[ Then we show that a unique generalized quasi- 
momentum satisfying properties 1-7 in Lemma 13.101 exists for general 
quasiperiodic solutions. Finally we check, that the motion is periodic 
if and only if the properties 8-9 are fulfilled. 

Assume, that the trajectory is periodic with the period T. Denote 
the impact points by ti,. . . ,tM , < ti < ... < Im < T. Consider the 
following function: 

nn N_ ^(7,^1 -0)^(7,^2-0) ^(7,^) 

^ ^ "^^^ ^(7,0) *(7,ti + 0)">(7,tM + 0) 

Lemma 3.7. The function 2(7) has the following analytic properties: 

(1) 2(7) is holomorphic in F outside the points A = and impact 
points Df. 

(2) 

iT 

(32) S(7) = 1 + — + ..., as X^oo 

V A 

(3) There exists a constant L such, that 

is regular near A = 0. 

(4) 2(7) has a pole of order at the point D°^^ and a zero of order 

at the point where denotes the number of impacts to 
the reflector Q HQs- 

Let us introduce generalized quasimomentum p for the billiard mo- 
tion. We would like to stress, that the generalized quasimomentum 



(33) exp(-4^S(7) 
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does not coincide with the "true" quasimomentum for the non- 
smooth function u{t) describing the billiard motion. The "true" quasi- 
momentum is defined on an infinite-genus Riemann surface. 
By definition: 



The function ^(7) is multivalued in F but its differential dp is well- 
defined. 

Lemma 3.8. The function ^(7) is a single-valued function on T* = 



r\ U % )■ ^''^ other words T* is obtained from T by removing the 



divisor trajectories. 

The proof of the Lemma follows immediately from the formula: 



The formula is well-defined if the denominator does not vanish. 

Lemma 3.9. Assume, that the gap [62^-1, &2i] contains at least one 
reflector, and that 7 G F lies inside the gap but outside the divisor 
trajectory. Then Re(S(7)) = 0. 

The proof follows immediately from the formula (l35l) . 
Now we are ready to formulate the properties of the generalized 
quasimomentum. 

Lemma 3.10. (1) dp is a meromorphic differential in T. 

(2) dp is holomorphic outside the points A = 0, D]^, 

(3) dp has first-order poles at the points D™, /}out ^^^^ 

(36) Res dp = — Res dp 

J^in £)out 

If we have two reflectors ds, ds+i in the same gap, then 

(37) Res dp = Res dp. 

nin r)in 

(4) dp has a second- order pole with zero residue at \ = 0. 



(34) 



P = ^ In (2(7)) 




(35) 




(5) 



(38) 




as X 



00. 
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(6) 

(39) ^dp = 0, j = l,...,n. 

(7) 
(40) 

^ dp = 0, if the gap [b2j-i, b2j] contains at least one reflector. 
Pi 

Let us recall that for gaps with reflectors the cycles j3j do not 
intersect the divisor trajectories. 



(41) Resdp=^, m, e N 

Di° iT 



(9) 
(42) 



where m.s denotes the number of impacts to the obstacle Q\^Qs 
forO<t <T. 



/27m ■ 
dp = ^ , Uj G Z if the gap [62^-1; &2j] contains no reflectors. 

Properties [6|, [7] follow immediately from Lemmata I3.8|3.9l and the 
fact that these cycles do not intersect the divisor trajectories is essen- 
tial. All other properties are completely standard. 

Remark 3.11. Differential dp is even with respect to the transposition 
of the sheets ofV. Therefore the second term in ^3^) can be estimated 
more accurately: 

(43) dp = d (y^^=j + O {^^J^ d\, as A 00. 
At X = one has the following expansion: 

(44) dp=^d(^-^^+Oil)d{V\), 
where L is the same as in /E 



Let us have a generic spectral curve T with a set of reflectors dg (with 
the reality conditions, formulated above). 

Lemma 3.12. There exists an unique differential with the properties 

UT 
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This differential plays the role of the generalized quasimomentum 
differential for generic billiard spectral data, and we shall denote it 
with the same symbol dp. 

Proof of Lemma 13.121 Let us construct dp as the following linear 
combination: 

(45) dp = Cdp' , where dp' = dp + ViSlfj^-^, 

J denotes the set of all a-cycles, containing reflectors, dp is the Abelian 
differential of the second kind with a second-order pole at the point 
A = such, that 

(46) dp = d{—=\^... as A — s> 0, (f dp = 0, j = l,...,n, 



are the Abelian differentials of the third kind with first-order poles 
at the reflector points in the a^, 

(47) Res fi(fe) = +1, Res ^(fc) = -1, \ D^""*) G a^. 



(48) ^fi(fc) = 0, j = l,...,n. 



From the Riemann bilinear relations it follows, that 

(49) Bik = j^i^k) = -2m j ui. 

By Corollary 13.51 the matrix Bik, /, A; G is non-degenerate, therefore 
conditions fHOl) defines the constants Vk uniquely. 

Let us consider the divisor dynamics with reflections in F, written 
in terms of the KdV time a (see (fT4l) ). Then outside the divisor tra- 
jectories we have: 

L 

(50) dp' = Urn — [ ^^da 



For the logarithmic derivative of the wave function we have the follow- 
ing formula 



(51) 



2n 

^'(7, a) p=i w + \Rn-iiX) 



^(7,^) f^x.ia) A(A-Ai(a))...(A-A„(a))' 
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where i?„_i(A) denotes some polynomial of degree n — 1. For A — >^ oo 
we have 



(52) 




i=i 

therefore the leading term at oo does not vanish and has a constant 
sign. As a corollary the leading term of dp' at oo does not vanish, and 
the constant C is defined by fl38|) . 
Lemma 13.121 is proved. 

Let us formulate the sufficient conditions for the periodicity of the 
billiard trajectory. 

Proposition 3.13. Assume that the quasimomentum differential dp, 
defined by the properties\T^ satisfy also properties\^ [5|/or some T . 

Then the corresponding divisor billiard motion is periodic and with 
the period T in the natural parameter, and there exists an integer K 
such, that the corresponding billiard trajectory on the ellipsoid is peri- 
odic with the period KT . 

Let us consider the divisor dynamics with reflections in F, written 
in terms of the KdV time a (see (1141) ). Denote by A{t) the Abel 
transform of the divisor expanded by the auxiliary basis cDi,. . . , u;„. 
From a standard calculation based on the Riemann bilinear identities 
it follows that 



dA {a) \ ~^ if the gap [62j_i, &2j] contains reflector points 
^ if the gap [62j_i, &2j] contains no reflectors. 



da 

where L is the same as in (l33l) . (I44l) . 

Let us denote by Qi the union of all trajectories generated by the 
divisor motion for < cr < L and by Qq the following collection of 
oriented curves: 

u "a|u( u 

, all gaps containing reflectors ) \ all gaps without reflectors 

From the Dubrovin equations (j5]) it follows that the direction of motion 
in each gap is constant. Therefore the difference between Qi and 
satisfies the conditions of Lemma 13.31 The Abel transform for this 
difference is zero, therefore these two sets of paths coincide. It implies 
the periodicity of the divisor motion with the period L in the variable a. 



r 
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Following the procedure described above one can naturally introduce an 
analog of the generalized quasimomentum with respect to the variable 
a. Let us denote it by rfp'. It coincides with d'p up to a normalization 
constant, and it is easy to check, that lip' = ^(ip, therefore the divisor 
motion is also periodic in t with the period T. 

Remark 3.14. Let us show, that both of the periodicity conditions 
formulated above are equivalent. Equation (d^j defines a collection of 
numbers mj uniquely up to a constant multiplier. From the bilinear re- 
lations we see that they are proportional to either residues or f3-periods 
of the generalized quasimomentum differential dp whether or not the 
corresponding gap contains reflectors. Proposition \3.13\ means exactly 
that in the periodic case one can choose the constant multiplier so that 
all rrij are integer numbers. 

4. ISOPERIODIC DEFORMATIONS FOR THE BILLIARD MOTION 

Denote the zeroes of the quasimomentum differential dphj rji,. . . ,rin+r- 
Then the quasimomentum differential takes the form: 



(53) dp 



By analogy with [22] we shall use the following statement: 

Lemma 4.1. Assume, that we have a periodic spectral data for the 
billiard motions and a deformation of these data such, that 

d dp 




(54) 

or, equivalently, 
(55) 



dp 



A=const 



= ^(7), 

A=const 



where -P(7) is a meromorphic function in T . Then this deformation 
preserves the T -periodicity of the divisor motion associated to the bil- 
liard dynamics. 

In our situation the spectral curve depends on the variable C- While 
calculating the partial derivatives with respect to C, we assume the 
projection of the point 7 to the A-plane to remain constant. We remark 
that this rule generates singularities at the moving branch points. 
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Remark 4.2. The function p in the left-hand side of ( f53]) is multival- 
ued. The condition that its derivative is single-valued is equivalent to 
the isoperiodicity condition (see [22] ). 

The proof of the Lemma is rather straightforward. From the equation 
fl54p it follows that all periods of dp as well as the residues at the 
points D™, are constant, therefore the properties [HI E] are satisfied 
identically. 

Lemma 4.3. Assume, that we have a deformation of the periodic spec- 
tral data satisfying (54\ ). Then -P(7) has the following analytic proper- 
ties: 

(1) P{'~f) has 1-st order poles at the points A = and D™, 1^°"*. 

(2) P(7) has 1-st order poles at the finite ramification points hi,. . . ,h2n- 

(3) -P(7) has a third-order zero at X = oo. (Here we use that dp is 
even with respect to the hyperelliptic involution t). 

(4) P(7) is odd with respect to the hyperelliptic involution r: P{t'j) = 
-P(7). 

The proof again is rather straightforward. We prove now that the 
opposite is also true. 

Lemma 4.4. Assume that P{'j) is an arbitrary meromorphic function 
satisfying the properties [2H3 from Lemma \4.'^ Then it generates an 
infinitesimal isoperiodic deformation. 

This lemma generalizes the analogous statement from [22]. The idea 
how to construct this deformation is very similar to the arguments of 
the paper [20] . 

First of all let us define the dynamics of the branch points. Consider 
the expansions of p, P near the point hj-. 
(56) 



P = Po,b,iO+Pi,b,iO^^- bj{0 + - ■ ■ dp 



2VA-&,(C) 



dX 



(57) P = ^i^ + ... 



1 ^ rn^ + 



Comparing (l57|) with ( !58l) we obtain: 

•^C Pi,bAC) 
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Let US define now the dynamics of tlie reflectors dg. Near tlie points 
D™, D"'^^ we have respectively: 



(60) p = ±-^ln(A-4(C)) 



dp 



-zT{X-d,{0) 



dX + ... 



(61) 



(62) 



P 



(A -4(0) 



dp 



dd.io 



A=const ^T(A-4(C)) 

Comparing (16T1) with (!62l) we obtain: 



+ 



(63) 



Consider the deformation defined by fISU]) . fl55]) . Then equations 
f lSS]) define a deformation of ^(7) as well as a deformation of rfp such, 
that the deformed function dp is defined on the new spectral curve and 
satisfies all properties of Lemma I3.10[ Therefore the deformed spectral 
data also generate periodic divisor motions with the same number of 
impacts. 

By analogy with [22] it is also convenient to derive the dynamics for 
the zeroes of dp. In the neighborhood of rij we have 



(64) 
(65) 



# = ±Pi,^,(C)(A-r/,(C))rfA + ... 
P = ± [Po,.,(C) + i^i,.,(C)(A - vAC))] + • • • , dP = ±Pi,r,,iC)dX + ... 



(66) 



(67) 



d dp 



A=const 



Pl,V, (0 



dX + 



Vl,ru (C) 



We have shown that we have a one-to-one correspondence between the 
period preserving deformations and meromorphic functions with special 
analytic properties. By analogy with [22j let as introduce a convenient 
basis in the space of such meromorphic functions. Let us define: 



(6J 



(A — rij)dX' 



3 



,n + r. 
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In the paper by Marchenko and Ostrovski [31] the following coordinates 
were used on the space of spectral curves associated with periodic po- 
tentials: 

(69) (j = -ip(nj), j = l,...,g, 

where r]j are the zeroes of the quasimomentum differential. It is conve- 
nient to use these coordinates in the method of isoperiodic deformation. 
In particular, the corresponding flows commute at regular points for 
these flows (all zeroes r]j are pairwise distinct). 

This choice of local coordinates correspond to the following choice of 
the factors Cj: 

(70) PAVJ)=^■ 
We use here, that 

dPiVj) 



A=const,7=r7j 

at the zeroes of the generalized quasimomentum differential dp. There- 
fore we obtain 

2n r 

-Vj n ivj - bk) n ivj - dsf 

It is convenient to assume Cj > for real zeros of quasimomentum 
differential. 

Calculating the leading terms of Pji'y) at the branch points, reflectors 
and zeroes different from rij we obtain: 

dbk ^ bk 



dQ 
dds 
dQ 



-'3 



■3 



bk 


-Vj 




ds 


ds 


-Vj 




Vk 


Vk 


-Vj 



(7# = C,—^, k^j 



^Vj r 
8Q - '^^^ 



L_l V 

r,. - 71, 9t) . 9 2-^ r,. -h, 2-^ 



rij-rik 2rij '^f^^Vj-bk jr^Vj-ds 



-k^j 

Summarizing all these results we obtain: 

Theorem 4.5. Let us consider a generic set of spectral data corre- 
sponding to a periodic billiard motion (generic means that all zeroes of 
dp are pairwise distinct). Then equations (7^ with the normalization 
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locally define a set ofn + r commuting flows on the spectral data, 
preserving the T -periodicity of the billiard motion as well as the number 
of impacts during a period. 

Remark 4.6. In the billiard motion it is critical to respect the reality 
conditions. To keep the reality of the spectral data it is sufficient to 
assume, that the times (j corresponding to real zeroes rjj of dp are real, 
the times corresponding to complex conjugate pairs of zeroes rjk, rji =fik 
are complex conjugate: = Cfc- 



5. Appendix. Some details of numeric check. 

To illustrate the approach developed above, some numerical simula- 
tion were made. For computational purposes it is more convenient to 
use the spectral parameter E = 1/X instead of A. Equations fl73|) were 
implemented using the standard 5-th order Runge-Kutta method. The 
authors used the integrator, developed by Hairer, N0rsett and Wanner 
|24j . this free software package in available on the Ernst Hairer's web 
site http : / /www . unige . ch/~hairer/. The program was developed as 
a modification of the program used in [22j to check isoperiodic approach 
to the KdV equation. 

The novelty of the problem was, that the calculation of the starting 
point for the billiard problem required some additional calculations. 
Let us discuss the search of the starting point. 

It is convenient to calculate the generalized quasimomentum differen- 
tial for a small periodic perturbation of the constant potential u{x) = 1. 
We assume the period T = 1. In our experiments we assumed that we 
have one reflector in the last gap [&2n-i, &2n]- In our program the branch 
points are enumerate in the inverse order Ei = l/62n,- ■ • , E2n = 
For numerical reasons it is critical to make the gap, containing the re- 
flector, a little bigger than the other ones. We also assume, that at the 
starting point the reflection wall lies in the center of the distinguished 
gap. The starting lengths are: E2 — E1 ~ 10~^ and E2k — E2k-i ~ 10~^, 
k>l. 

First of all, let us calculate the positions of the point r (the center 
of the gap), and the zeroes Aq = l/rjo, Ai = l/r/i (up to some higher 
order corrections). 

Let 
(74) 

Ei = C-A, E2 = C+A, l/rf = R = C, Ao = C+Ao, Ai = C+Ai, A<1. 



PERIODIC BILLIARD ORBITS ON n-DIMENSIONAL ELLIPSOIDS ... 

Assuming all other gaps to have zero length (genus 1 approximat 
we have 

duji — H — ^=dE 
y/E{E-Ei){E-E2) 

Let us denote 

E^Z + C 

The constant x is defined from the condition 

Z=A 



z=o 

Expanding dcui near Z = we obtain: 

1 dZ X / Z ,A dZ 

dui = X , — = ^ 1 + 0(Z^) , 

Using the table integrals 

Z=A 

/dZ 7ii 
VZ?^ " ~2 

z=o 



Z=A 

ZdZ 



we obtain: 



2=0 

Z=A 



VZ2 - A2 



iA 



Therefore 



z=o 

- c) 

In the genus 1 approximation we have 

2 - R)^yE{E - Ei){E - E2) 
Let us look at the conditions on Aq,Ai imposed by 

dp^O 
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Expanding near Z = we obtain 

'(Z-io)(Z-ii) \ dz 



z \ vz2 - A2 



B ( A^Ax /r T AqA^ ^, \ dz 
2v^ ^ Z y 2C J ^ ' VZ2 - A2 



Due to parity constraints 



7 Z VZ2 - A2 7 VZ2 - A2 



and finally we obtain: 

(75) i, + A, + ^ = 0(A^) 

It is difficult to use the condition 

dp = Q 



directly. It is more convenient to use Riemann bilinear relations in- 
stead. We know, that if 



then 



dp = <p dp = 0, 



Res dp = RespdJi 

Rin OO 



Due to the periodicity condition 



I.e. 



But 



Res dp = —i, 



Resp Ui — —i 



therefore 



(76) Bx = I 



B = ^ ^ = ^^+0(A^) = ^^+0(A^ 

C'' + C{A + A^)+AoAi C2 + ^ 1 + ^ 
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Direct calculation of the residues of dp at the points i?™,/?"*^*, gives 

us 

1 = 

2^/Ck 

Therefore 



(77) ioii ~ -2aVc 1 + ~ -2Av^ 1 



\^ 2C2 y V Cv^ 



To determine C let us use (1761). 



1 1 

X 



1 1 

X 



Finally, we obtain 



(79) C 



To find the starting point we first solve (1791) . then we calculate Aq, 
ii using (I77D, (I7HD 

Let us calculate the resonant points, defined by the condition 

p(Cfc) = nk, k > 1. 

To do it let us expand dp outside the neighborhood of the point C. 

B dZ (Z - y4o)(Z - Ai) 



dp , , 

2 v/zTc - A2 

B dZ if AqAi 



2 VzTc - A2 y z 



(io + v4i) + Z 



2 yzTcVZ (^-(-^oH-^OH-Z 

B ^Z / _ ip + ii ioA , ^.,2 
2v/Z + cl^ Z + Z2 

-^^^l-^^^ + ^^ + OfA^) 

~2^r E-c ^{E-cy^^^^> 
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Using the table integrals: 



dE 



dE 



1 { Ve-Vc 

log 



2VE{E-C) 2v^ \VE + Vc 



dE 



2y/E{E-Cy 4C 



2y/E 1 fy^E-Wc 

E-c^ vc^'^yvE+vc 



we get 



p = B 



E-\Ao + A,+ 



AqA^ 



2C I 2VC 
Taking into account fl75l) we obtain 



1 , / y/E-y/c\ AqAi y/E 
log 



E + VC 2C E 



p = bVe 



2C E-C 



Therefore we calculate the resonant points from the equation 

I 2 

irk 1 

^0) 



B 1 



2C(Cfc-C) 
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